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HYDRODYNAMIC  ASPECTS  OF 
THE  SPLIT-FOIL  LASER  TARGET  DESIGN 

I.  INTRODUCTION 

The  generation  of  uniform  electron  number  densities  in  neon-like  states  of  high-z  plasmas 
which  develop  from  exploding  foils  has  become  a  topic  of  much  recent  interest.  Such  plasmas 
are  necessary  for  the  creation  of  a  pumped  soft  x-ray  laser.  To  date,  perhaps  the  most  successful 
experimental  laboratory  x-ray  laser  has  been  the  one  designed  and  operated  by  Matthews 
and  co-workers  at  Lawrence  Livermore  National  Laboratory  (Matthews  et  al  1984).  In  these 
experiments,  an  amplified  x-ray  beam  is  formed  when  an  selenium  foil,  which  is  irradiated 
from  both  sides  for  0.6ns  by  0.53 /xm  laser  light,  explodes.  A  uniform  cylindrical  plasma  with 
electron  number  densities  and  temperatures  appropriate  for  x-ray  lasing  is  thereby  generated. 

It  is  of  interest  to  generalize  these  results  to  other  systems  and  parameter  regimes.  Of 
particular  promise  are  the  conditions  obtainable  with  the  Pharos  m  laser  at  the  Naval  Research 
Laboratory.  In  one  mode  of  operation,  this  laser  can  provide  a  single,  500  Joule,  1.06/im  beam, 
with  rise  time  on  the  order  of  1  ns,  and  peak  laser  intensity  of  approximately  1  x  1013~£r- 
Considerable  investigation  into  the  use  of  this  system  to  generate  a  pumped  soft  x-ray  laser 
has  been  undertaken  (Elton  et  al  1986;  Lee  et  al  1986).  We  extend  this  research  by  means  of 
numerical  simulation,  with  initial  conditions  and  parameters  appropriate  to  the  NRL  Pharos 
III  laser  operating  in  its  single-sided  illumination  mode. 

We  present  results  from  the  exploration  of  a  somewhat  idealized  NRL-like  set  of  regimes 
in  which  a  measurable  x-ray  beam  can  be  formed.  Sec.  II  is  a  discussion  of  the  geometry 
considered,  the  assumptions,  and  the  numerical  methods  used.  In  Sec.  HI  we  discuss  re¬ 
sults  from  the  simulation  of  the  hydrodynamic  evolution  of  solid  copper  foils  irradiated  by  a 
1.06/xm  laser  beam,  [CASES  A  -  Fj.  Results  from  these  simulations  indicate  that  the  evolved 
electron  number  densities  decay  from  an  upper  bound  of  approximately  3  x  1020  particles/cc, 
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at  approximately  0.25  ns  after  burnthrough.  A  promising  alternate  foil  design,  the  split-foil, 
is  then  presented,  in  Sec.  III.  It  is  found  that  relatively  uniform  plasmas  with  electron  number 
densities  of  6  x  1020  particles  per  cc  can  be  generated  with  1/xm  thick  split  copper  foils,  with 
a  slit  on  axis  of  width  of  0.01  cm.  A  discussion  of  these  results  is  then  given  in  Sec.  IV. 

II.  SIMULATION  MODEL 

The  described  simulations  are  performed  with  the  FAST2D  laser-shell  simulation  code 
(Boris,  1976;  Boris,  1977;  Emery  et  al  1982).  The  code  solves  the  equations  of  ideal  compress¬ 
ible  hydrodynamics  in  two  full  dimensions,  using  the  flux-corrected  transport  algorithm  (Boris 
and  Book,  1976).  In  this  code,  energy  is  deposited  in  the  plasma  by  the  incident  laser  beam 
through  classical  inverse  bremsstrahlung  absorption,  and  transported  through  the  plasma  by 
classical  nonlinear  thermal  conduction.  An  adaptive  zoning  algorithm  is  employed  to  insure 
sufficient  resolution  in  the  regions  of  strong  density  gradients. 

The  geometry  is  identical  for  all  the  simulations  here  described;  see  Fig.  1.  The  laser  target 
is  assumed  to  be  much  longer  in  the  direction  of  eventual  x-ray  lasing,  the  z-direction,  and 
uniform  in  that  direction.  This  allows  the  computation  to  be  limited  to  two-dimensional  cross- 
sections  of  the  exploding  foil  and  to  be  accurate  without  sacrificing  efficiency.  In  addition, 
the  problem  is  assumed  to  be  symmetric  about  the  “vertical”  axis,  y.  The  equations  thus 
only  need  be  time-advanced  in  the  upper  y  half-plane,  (the  shaded  region  of  Fig.  la,  b), 
with  perfectly  reflecting  boundary  conditions  imposed  at  y  =  0.  Boundary  conditions  which 
simulate  unimpeded  outward  flow  of  the  plasma  are  imposed  at  ymax ■  Outflow  conditions 
also  are  imposed  at  xmax,  and  at  xmjn.  By  means  of  an  adaptive  gridding  algorithm,  large 
buffer  zones  are  constructed  at  both  xmax  and  xm<n;  these  buffer  zones  isolate  the  important, 
central  region  of  the  mesh,  and  effectively  set  xe<ige  to  ±oo.  The  laser  beam  is  assumed  to 
originate  at  x  =  -foo  and  the  target  is  centered  at  x  =  0.  Typical  runs  are  performed  on  grids 
with  120  points  parallel  to  the  incident  beam  (x-direction),  and  40  points  in  the  transverse 
direction  (y-direction).  In  order  to  resolve  the  large  density  gradients,  many  fine  zones  of  the 
non-uniform  mesh  are  clustered  in  the  neighborhood  of  the  regions  of  large  density  gradients. 
Since  a  variable,  Courant-limited  timestep  is  used  in  the  calculation,  the  timestep  is  thus 
generally  small,  of  the  order  of  0.02  ps.  The  code  takes  approximately  0.1  second  per  timestep 
of  CRAY-XMP  time,  or  about  15  minutes  for  each  simulation  considered  here. 


Initial  conditions  for  both  the  solid  foil  simulations  and  for  the  split-foil  runs  are  specified 
by  the  one-dimensional  quasi-static  equilibrium  model  (Orens,  1980;  Manheimer  et  al  1982). 
In  the  split-foil  cases,  the  gap  is  similarly  modelled,  for  code  stability,  with  initial  density  (x 
=  0,  y  =  0)  less  than  critical. 

In  each  case,  the  material  chosen  for  the  foil  is  pre-ionized  to  the  neon-like  state.  From 
evolved  densities  and  temperatures  we  expect  that  this  is  a  state  the  plasma  would  attain 
(Rosen  et  al  1984).  In  addition,  the  ionization  level  is  held  fixed  throughout  the  simulation. 
Ijiis  constant  Z  assumption  is  computationally  efficient.  Since  the  timescales  for  ionization 
and  recombination  vary  as  the  inverse  of  the  density,  becoming  much  longer  than  relevant 
hydrodynamic  timescales  in  the  expansion  phase  of  the  generated  plasma  (London  and  Rosen 
1986),  the  assumption  is  also  physically  appropriate. 

We  now  turn  to  a  discussion  of  individual  foil  cases  and  results. 

III.  SOLID  FOIL  RESULTS 

Table  1  lists  the  runs  we  have  chosen  to  display  here  to  demonstrate  density  and  temper¬ 
ature  states  that  have  evolved  from  a  range  of  thin,  solid-  copper  foils. 

For  all  simulations  of  this  section,  the  incident  laser  wavelength  is  1.06/xm,  with  a  peak 
intensity  that  is  linearly  ramped  up  for  1  ns,  followed  by  a  flat-topped  peak  intensity,  specified 
for  each  case.  Spatially  the  pulse  is  Gaussian,  with  full-width  half-maximum  100/xm  from  the 
center  of  the  foil. 

Since  the  time  of  0.25  ns  after  burn-through  has  generally  been  found  to  be  a  time  by 
which  the  exploded  foil’s  sharp  density  gradients  are  no  longer  observable,  it  is  the  time  chosen 
for  comparing  CASES  A-F,  as  indicated  in  Figs.  2-7.  Although  a  desirably  uniform  electron 
density  profile  has  been  obtained  in  each  of  the  six  cases  by  this  time,  the  electron  number 
density  in  each  case’s  physical  plateau  region  has  monotonically  decreased  to  a  relatively 
low  value  of  3  X  1020  particles  /  cc.  Variation  in  the  initial  foil  thickness  alters  the  time  of 
burnthrough,  and  alters,  somewhat,  the  electron  temperature  in  the  uniform  density  region; 
however,  the  electron  number  density  evolves  to  a  state  which  is  nearly  identical  in  all  cases. 
See  Figs.  3,  5  and  6,  and  Figs.  2  and  7.  Keeping  the  foil  thickness  constant  and  varying 
the  laser  intensity  again  alters  burnthrough  time  and  temperature,  but  does  not  change  the 
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estimated  number  density;  compare  Figs.  2,  3  and  4,  and  Figs.  6  and  7.  We  conclude  that  for 
thin  copper  foils  that  will  burn  through  by  approximately  2  ns,  with  single-sided  illumination, 
the  resulting  electron  number  densities  are  approximately  3  x  1020.  This  result  is  obtained 
when  both  peak  intensities  and  initial  foil  thicknesses  are  varied.  A  plasma  with  a  higher 
plateau  electron  number  density  is  desirable  for  greater  x-ray  gain  (Feldman  et  al  1984). 


IV.  SPLIT-FOIL  RESULTS 

In  order  to  generate  a  pumped  soft  x-ray  laser,  two  qualities  are  required  of  the  electron 
number  density:  (1),  the  electron  plasma  must  be  of  sufficient  density  for  x-ray  lasing,  and  (2), 
a  region  of  uniform  density  must  exist,  to  minimize  refraction  losses  (Chirkov  1984).  Alternate 
foil  designs  are  thus  considered,  in  an  attempt  to  both  raise  and  enlarge  the  plateau  region 
electron  number  density.  The  simplest  case  that  can  be  engineered  is  a  thicker  foil  than  that 
considered  in  Sec.  II.  While  thicker  solid  foils  would  not  be  feasible  if  burnthrough  by  2  ns  is 
desired,  a  design  geometry  that  increases  the  plasma  density  above  that  of  a  flat  foil  that  is 
only  1500A  thick  is  required.  By  initializing  the  problem  with  a  1/xm  thick  foil  which  has  a  slit 
at  the  y  =  0  axis,  we  find  that  regions  of  electron  number  densities  that  are  both  uniform  and 
approximately  6  X  102°  particles/cc  can  be  generated.  In  addition,  the  plasma  which  results 
from  this  novel  foil  design  is  found  to  survive  in  near  steady-state  for  over  1  ns. 

For  the  split-foil  simulations  considered  here,  the  copper  is  pre-ionized  to  the  neon-like 
state,  exactly  as  it  was  for  the  solid  foils.  In  addition,  the  incident  laser  pulse  is  identical  to  that 
used  for  the  solid  copper  foil  CASE  E,  with  a  peak  incident  laser  intensity  of  1.00  x  1013~~y- 

Fig.  8,  plots  of  peak  axial  electron  number  density  measured  at  the  center  of  the  slit, 
ne(x,y  =  0),  as  a  function  of  time,  depict  the  essential  result  of  the  split-foil  simulations.  Fig. 
8a,  a  plot  of  peak  ne  (x,  y  —  0)  vs  time  for  a  split  foil  that  was  initially  10/zm  thick,  shows  a  peak 
electron  number  density  that  rises  above  critical,  to  stabilize  near  a  value  of  approximately 
1.3  x  1021  particles/cc.  Fig.  Sb  is  a  similar  plot  for  the  1/rm  split  foil.  In  this  figure,  it  is  clear 
that  the  peak  ne{x,y  =  0)  hovers  near  6  x  1020  particles/cc  for  longer  than  1  ns. 

Fig.  9a,  a  contour  plot  of  ne  from  the  split-foil  case  at  t  =  1.9  x  10~9s,  demonstrates 
how  ne  has  plateaued  in  the  initially  split  region,  with  approximately  6  x  10'°  particles  / 
cc.  The  plateau  region  is  approximately  100/rm  in  the  x  direction  by  50/xm  in  y  (40/jm  in 
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the  half- plane).  A  contour  plot  of  electron  temperature,  Fig.  9b,  shows  that  at  this  time 
the  temperature  in  the  split  region  is  approximately  600  eV.  Note  the  uniform  density  region 
across  which  a  uniform,  controlled  temperature  gradient  appears.  This  should  be  useful  for 
calibrating  the  lasing  models. 

To  check  the  sensitivity  of  our  results  to  material  variation,  we  also  performed  simulations 
with  selenium.  For  a  given  temperature,  the  ratio  of  the  hydrodynamic  pressure  to  the  mass 
density  essentially  depends  on  the  ratio  of  a  target’s  ionization  level,  Z,  to  its  atomic  mass, 
A.  Thus,  once  a  plasma  has  been  generated  and  heated  by  the  incident  laser  beam,  the  hy¬ 
drodynamic  evolution  that  plasma  undergoes  will  be  governed  by  the  ratio  Z/A.  For  a  rather 
wide  range  of  materials  of  interest  this  ratio  is  nearly  constant.  For  instance,  the  ratio  Z/A 
for  copper  ionized  to  the  neon-like  state  is  0.299.  The  ratio  Z/A  for  selenium  ionized  to  the 
neon-like  state  is  0.303.  From  this  similarity,  it  may  be  expected  that  results  from  simulations 
with  selenium  ionized  to  the  neon-like  state  will  be  nearly  identical  to  results  from  the  above 
described  simulations  with  copper.  We  find  this  to  be  the  case;  results  from  a  simulation 
performed  with  a  pre-ionized  split  selenium  foil  that  initially  was  1/zm  thick  are  very  similar 
to  those  obtained  with  copper  foils.  The  incident  laser  pulse  used  in  the  selenium  simulation  is 
spatially  uniform,  and  is  temporally  identically  to  that  specified  for  CASE  E.  Fig.  10a,  a  plot 
of  peak  axial  ne(x,y  =  0)  vs  time  for  the  split  selenium  foil  demonstrates  that  the  hydrody¬ 
namic  evolution  of  the  resulting  plasma  is  very  similar  to  that  generated  from  copper.  Fig.  10b 
is  a  snapshot,  at  t  —  1.125  X  10-9s,  of  the  electron  number  density  which  has  developed  from 
the  1/xm  thick  split  selenium  foil.  This  figure  shows  the  large  axial  region  of  uniform  density 
which  has  been  formed  from  the  ablating  selenium  foil,  a  result  which  is  nearly  identical  to 
results  obtained  from  copper  foils. 


V.  DISCUSSION 

We  have  performed  numerical  simulations  of  two  distinct  copper  foil  target  designs.  In 
all  the  cases  simulated,  the  exploding  thin  solid  copper  foils  with  single-sided  illumination 
generate  somewhat  uniform  electron  number  densities  which  are  approximately  bounded  from 
above  by  3  x  102°  particles/cc.  This  low  number  density  in  the  physical  plateau  region  is 
obtained  even  when  the  initial  target  foil  is  thickened  by  a  factor  of  two. 
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Plateau-region  plasma  electron  number  densities  can  be  varied,  however,  as  simulations 
with  ablating  split  copper  foils  demonstrate.  For  10/xm  split  copper  foils,  super-critical  axial 
densities  are  achieved.  From  ablating  split  copper  foils  that  are  initially  Ifx m  thick,  approxi¬ 
mately  100/xm  plasmas  with  relatively  flat  electron  number  densities  of  approximately  6  x  1020 
particles  /  cc  are  generated.  This  higher  density  plasma  is  found  to  exist,  in  near  steady-state, 
for  over  1  ns.  A  key  to  understanding  the  reason  for  the  development  of  a  nearly  steady-state 
situation  can  be  found  in  the  behavior  of  the  plasma  flow. 

Two  evolved  velocity  fields  are  depicted  in  Fig.  11.  Fig.  11a  is  the  flow  pattern  for  the 
solid  copper  foil  CASE  E  at  t=  2.12  x  10-9s.  Fig.  lib  is  the  flow  pattern  at  t=  1.9  x  10~9s 
for  the  lfim  split  copper  foil  simulation.  The  incident  laser  pulse  is  identical  in  each  case.  A 
dynamical  explanation  for  the  difference  in  peak  axial  electron  number  densities  which  results 
from  the  foil  shape  may  be  obtained  from  this  pair  of  figures.  In  the  solid  foil  case,  a  flow  field 
is  set  up  that  can  only  convect  the  plasma  rapidly  away  from  the  region  of  interest.  On  the 
other  hand,  from  Fig.  lib  it  is  clear  that  a  nearly  steady-state  situation  has  evolved  in  the 
split-foil  case.  As  on-axis  plasma  is  swept  away  to  x  =  ±oo  a  nearly  equal  amount  is  steadily 
ablated  into  the  neighborhood  of  the  axis  from  the  relatively  infinite  solid  copper  reservoir, 
the  split-foil.  The  pressure  and  density  are  higher  in  the  split  foil  case  because  the  vertical 
velocity  must  be  stagnated  and  because  the  mass  available  will  only  depend  weakly  on  the 
pulse  intensity.  This  steady-state  situation  is  found  to  be  robust,  as  is  shown  by  Fig.  8b. 
In  this  steady  situation,  the  plasma  which  streams  away  from  the  potential  lasing  region  is 
continually  replaced  by  plasma  that  is  fed  into  the  region  from  the  ablating  foil.  This  evolution 
is  also  observed  in  hydrodynamic  simulations  of  split  selenium  foils. 

Thus,  by  means  of  this  new  target  foil  design,  we  find  that  generated  plasma  densities 
may  be  controlled,  both  as  a  function  of  time  and  of  space. 
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<  a  >  Solid-foil  computational  domain 


Fig.  1  —  Simulation  geometry 


(b)  Split-foil  computational  domain 
Fig.  1  (Continued)  —  Simulation  geometry 
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(b)  Contours  of  constant  temperature,  in  eV 
Fig.  4  —  Thin  solid  foil,  CASE  C,  at  t  ™  1.46  x  10~9s 
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(b)  Contours  of  constant  temperature,  in  eV 
Fig.  7  —  Thin  solid  foil,  CASE  F,  at  t  —  1.8  x  I0“9s 
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(b)  Contours  of  constant  temperature,  in  eV 
Fig.  9  —  lum  split  copper  foil,  at  a  time  1.9  x  10-9s 
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(a)  Peak  axial  electron  number  density  as  a  function  of  time 
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(b)  Contours  of  constant  ne/10:o,  at  a  time  1.9  x  lCT’s 
Fig.  10  —  1/itn  split  selenium  foil 
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